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TECHNICAL NOTE NO. 797 



THE END-PLATE EFFECT OF A HORIZONTAL -TAIL SURFACE 
ON A VERTICAL -TAIL SURFACE 
By S. Katzoff and William Kutterperl 



SUMMARY 



The end-plate effect of the horizontal -tail surface 
on the vertical -tail surface has "been theoretically studied 
by derivation . of the flow corresponding to the condition of 
minimum induced drag and "by solution of the usual lifting-, 
line equations. The results obtained "or the two nethods 
were considerably different, indicating that the condition 
of minimum induced dra? (uniform normal induced velocity) 
is. not a close approximation to the actual condition. The 
derived span load distributions' are" given, together with 
equations and curves for the total lift forces and the mo- 
ments of the horizontal? and the vertical -tail surface's 
about the intersection. The results obtained" Indi cat e "a '" "" 
considerable increase in effectiveness .of the vertical ' "" 
tail. For example, the effective aspect ratio, in terms 
of the elliptical airfoil, is about 1.5 times the geometric 
aspect ratio for a tail ?roup in which the span of the 
horizontal tail is twice the span of the vertical tail. 



INTRODUCTION . . 

m 

*" " — _ 

0 -. ^_ 

The effectiveness of the vertical tail is known to 
be considerably increased by the presence at its base of 
the horizontal tail. In conjunction with the tail-'surFace 
studies in progress at the 1TACA, a theoretical study has 
been made of this end-plate effect for the idealized case' 
of an isolated tail group (fig. l). 

In the usual theoretical end-plate study (references; 
1 and 2), the problem is -simplified by. considering only 
the case of minimum induced drag; that is, the case of 
constant normal induced velocity along the span of the 
lifting surface with zero normal induced velocity along 
the span of the end plate. '3y this' restriction' (reference 
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3), the problem reducos' to that of the two-dimonsi onal 
flow past a figure which represents cf cross section of the 
wake infinitely far downstream; in tho present case, the 

figure Is I -shaped (fig. 2). 

As will he shown, however, the condition of minimum 
induced drag is physically impossible , for it requlre-s 
that tho end plate have infinite chord. Preliminary- stud- 
ios having indicated that the corresponding inaccuracy 
might he consi derablo when the mothod is applied to the 
typical tail group, a core general method in which tho 
usual lifting-line equations of the monoplane airfoil were 
applied was used to obtain a more accurate solution. This 
method is discussed in part I; — the method that as sumo s 
minimum induced drag is considered in part II. 



SYH30LS 

,, 1 - wako cross section 

H, D dimensions of _l_ obtainod by transformation 

h span of- vertical tail 

d sonispan of horizontal tail 

r ratio of vortical-tail st>an to horizontal-tail ppan 

(h/2d or H/2D) 



m 



a function of r< = ii+J^JS^Z**!) 



r 

S area 

c chord 

A aspect ratio 

I lift per unit length 

L total lift or cross-wind force 

P air den_sity 

0 Ti lift coefficient 
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M moment 

a ' angle of attack 

a slope of lift curve 1,}/ 4a o ) 

a 0 slope of section lift curve 

Y free-stream velocity - ■ 

vr induced, velocity 

Y circulation per unit length (spanwise) in wake 
p bound circula.tion _. .. . 

x, y coordinates of projection of tail group 

Xj^, y x variables of integration 

Complex variables: 

z plane of 1 (x, y) 

£ plans, of unit circle (r, 9) 

p auxiliary plane 

tf potential function ($ , \j>) ... .. ... 

Subscripts 

x point on X axis 

y point on Y axis •• • 

h ' vertical tail 

d horizontal tail 

e effective 

0 geometric 

1 induced 



— • 
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I. FINITE HORIZONTAL TAIL 
Method 



The an^Ie of attack of the verti cal -tail surface is 
assumed to he a 0 and that of the horizontal-tail sur- 
face (the end plate) is zero, the corresponding vortex 
system oeins as indicated in figure 3. The hound vortex . 
at the hase of the vortical" tail divides at the- int ersec- 
tion and the two halves continue, with, opposite rotations, 
alone; the horizontal tail. 

The normal~~induced velocity w x at a point P (0,y) 

of the horizontal tail (fie. S) iua to the trailing vortex 
system is, hy the Biot-Savart law, 

f ^ 1 r y «** A . 1 r & 1 ard at 

w x<y> ■ ^ J dxl dXl + "to J y7~^7 dy, ^ 

0 -d (1) 

Similarly, for the normal induced velocity w y et a point 
Q(x,0) of the vertical tail, 

1/ x d% i r i dr h 

w y (x) = / g ~ dy x + / ~ dx x 

y 4tt x " + y x a dy x 1 47T z » "* x dl i 

From the relation "between the hound circulation and the 
effective anfle of attack, two additional equations are 
o htained : 

r d (y) = ^ (y) (3) " 

r h (x) = [va 0 + w y (x) J (4) 

r 

Derivation of the span load di strihuti on for any given 
case involves the simultaneous solution of equations (l), 
(2), (3), and (4). The most practical method found for 
effecting the solution was a graphical method of succes- 
sive approximation. A distribution of and Th * s 
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assumed and the operations of equations (l) and (2) are 
performed for a number, .of values of "x and of y to find 
the normal induced ve'lo cities w x (y) and Wy(s) at these 

points. Substitution of these values in equations (S) and 
(4) leads to_ a span ■ load, distribution that, in general, 
differs from the assumed distribution. The difference be- 
tween the two di stributions. will indicate what 'changes' 
should be made in the assumed distribution £o make the 
corresponding • derived distribution agree more closely with 
it.' In the application of this method tie process las 
found to converge very rapidly, so that very little addi- 
tional computation was required after the first or. the 
second approximation to get accurate agreement". 

Results 



Span load distributions were calculated for tail 
groups' in which the horizontal and t'he vertical tails js^re . 
of elliptical and semielliptical plan form, respecti vel'y , 
and have the same maximum chord. The value of a 0 was 1 ' 
assumed to be 5.3 per radian. The results are shown ih^ 
figure 4 for various aspect ratios and" for various ratios" 
of vertical-tail to horizontal-tail 'span r. "_ — ' *' -'" : 

At the intersection of the two tails, the slope of": the 
span load curve is positive for the. verti cal tail and (neg- 
atively) infinite for the horizontal tail. .This arranger^; 
ment can be shown to be the only one consistent with the 
condition that the bound vortex at the base of the vertical 
tail divides into two equal parts and continues aLong the two 
halves of the horizontal tail with opposite rotations. 

The total lift forces on the vertical tails were ob- 
tained by integration of the curves of figure 4. From 
these integrals were derived values of the effective as- 
pect ratio- A e , that is, the actual aspect ratio of ah 
elliptical airfoil having the same lift-curve slope, ac- 
cording to the equation 

|jj = a * a Q (l 
o 

The results are summarized in figure 5, where the 
ratio A e /A is plotted against r~~for various' values of 
A.- The curve for the condition of minimum induceST cCragV 
which is derived in part II, is included in the figure. 



a. 



) 
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The moments of the horizontal- 'and the verti cal -tail 
loads about the point of intersection were also obtained 
from the span load distributions; they are plotted in fig- 
ure 6(a); The horizontal -tail moment contributes appre- 
ciably toward reducing the restiltant moment. The, result- 
ant rolling moments about -the intersection are ^iveh in 
figure 7(a). Por comparison the corresponding moment 
curves for the minimum ihduced-draf? condition, which are 
derived in part II, are ^iven*in figures 6(b) and 7(b). 



II. TAIL SURPAC3S OP MINIMUM INDUCED DRAG 



As stated in the introduction, the case of the tail 
havin? minimum induced dra^ is solved when the two- 
dimensional flow past a _J_ -shaped figure Cf-i*?. 2) is de- 
termined. The tangential -velocity distribution at such a 
figure determined the bound vortex system, of the tail, 
which in turn determines the lift distribution. 

The Plow Past a -Shap'ed Pi^ure 

Outline of method .- The potential function of the 
flow past a JL. -shaped figure can be derived from the 
known potential function of the flow past a unit circle 
by means of the function that transforms the space, around 
the circle (£-plane) to the spaoe around the {sv-plaao). 

The transforming function z = f(£l_ must satisfv the fol- 
lowing conditions: 

(a) As £ .traces the unit circle, z traces 

the J_ 

(b) Lim z- = £ 

| — » 00 

(c) Por all points outside the unit circle, the 
transformation is one to one and confoTmal. 

The function z => f(£) can be derived in three steps: 

1. By means of a S chwar z - Chri stof f el transforma- 
tion (references 4r and 5) the real-axis of an auxiliary 
p-plane is transformed • into the I ■ , the upper half 
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of the p-plans ^oin<? into the region around the I 
(the z-plane ) . 

P.. By the known, function v ' 

• ' P - i 



the real axis of the. p-plane is transformed into the 
unit circle, the upper half of the p-plane soin<? into 
the region outside the unit circle (the "£-plane ) . 



3. She elimination of p from these two func- 
tions ^ives the desired transformation. 

Tra nsfo rmation of .the p-plan e to the z-nlaae . - By 

the S dhwarz-Chri stof f el theorem, the real axis of the p- 
plane is transformed into a _1 "by the function whose 

derivative is 

|j* Ki(p-u x ) " (p-u s ) " . . .(p-A) S (p-A)~ S • (6) 

where jj. x ., p, a , . . . are the exterior angles of the : 1 
(fi?5» 8), measured positive counterclockwise; u lt u a , 
. . . are' the points on the- real axis of the p-plane 
that transform to the. vertices of the' — L_ ; K t is a con- 
stant; A*. i's_the point in the p-plane that «;oes into 
z = <» ; and A is the conjugate of A. The vertex of the 
L thst corresponds'-tb ' p- f » is not considered in ex- 
pression (o). In accordance with equation (5) and condi- 
tion (b), it is clear that A = i. Also, "by symmetry, it 
is permissible to choose ~u x =.u 4 = a, -u 2 = u 3 = b, 
and u e = <»» Correspondini; points in the z- and the p- 

planes are shown in figures 8 and 9. Equation (6) thus 
becomes . • 

|» = g *Lz JL. (7) 

d ? y p a - a 3 ( P s + D fi 



which integrates to 
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pVp 3 - a 2 p/a 3 + 1 + y P a - a 3 
K l B — + I x 0 lo^ . — — ' 7 + 2 a 



p s + 1 p /T 3- ^! - 7? 3 - * 8 ~ 



1 1 + b 



(8) 



where B = — 

1 + a 



2 - • 3 



0 = 



■, a s . a o v 2 
1 a - a o - <£b 

4 ' (l + a a ) 3/3 



Equation (8) contains six eroitrary constants; a, b, and 
the real and tho imaginary parta of Ej. and K a « .--Cheso 
constants aro evaluated as follows: 

The coefficient 0 of the logarithm t-trrm.^«.st 
bo zero in order that z be single valued. By- this 
condi tion , 

b 3 = — ^-g . O) 

2 + a - 

The coefficient B now becomes a function only 
of a : 

3 = — 1 



2 + a 3 



The point of— intersection of the _L is placed at tho 
origin of the z-plane by putting z = 0 for p = i"** 
By this condition, 

K 3 = 0 

Equation (3) has , at this point, boon simplified to 
tho form 



z = 



a a 
2 + a p + 1 



< (10) 



A substitution for p from equation (5) ^ives 



Ki 7 1 + a a j + 1 / ,s „ /I - a 3 



i s + 2 (^Alt + i en)' 

XL + a ' 
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from which the value of K x is determined "by the condition 



that lim z = ti 



B 4 (2 + a s ) 
«/ 1 + 



The parameter K x is thus roal and positive, and the ori- 
entation of the JL is accordingly as shown in figure 10. 



fSivos 



A substitution for K x in equations (10) and (ll) 



z = 



1 + a' 



p -/p 2 - a 
p 2 + 1 



or, with 



m = 



l_ r _aj 
1 + a { 



(12) 



* * = i-i-i /V + 2 (i^|)t I 1 (18) 
£> v 1 + a 3 ^ 



(14) 



Finally, m oust "be determined so that the I 

will have the specified proportions. The length D (cor- 
responding to the semi span of .the horizontal tail) is 



found "oy substituting p = "fa = 



in equation (12). 



J 2 + a 

The length H (corresponding to the span of the vertical 

■> 

(a) 

H = 2-/2 (1 + m) \ (15) 

E " ' 



tail) is found by substituting p = oo: 

D = 1 - m 



r = 



2D 



y 2 (i 

1 - m 



(c) 
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From equation (15c), m may be found in terms of r: 



r s + l-*/l + 4r a 
m = 



Thus, as r varies from o to » (infinite end plate to 
no end plate), m rarios from -1 to 1. All the parameters 
of equation (8) have new been determi-rnrd . 

i G \ 
The unit circle . £ = e becomes, by equation (14), 



* g , 

e + 1 / 128 , „ i8 
TE — Ve + 2m e 



+ i 



e 



= 2 yiTcos I J cos 6 + m (loi 

which relates points on the unit circle to points on the 
— L . Corresponding points of the two figures are shown in 
figures 10 and 11, The stagnation points of the vertical 

flow past the circle ( Q & ± transform, by equation 

(16), to a - 2 >J m, corresponding to two co incident 
points on the real axis for m > 0 (r > «/ 2 ) and to two 
.symmetrically situated points ©n the imaginary axis for 
m $J 0 (r <</2~, the more usual case). 

T he -potential functio n o f the flow .-* The potential 
function of- the vertically upward flow past the circle 
with unit velocity is 

W = 1 V =.-1 (i - 

for points on the circle 

# = 2 sin 9 (17) 
\|/ = 0 

Equations (16) and (17) suffice to determine the velocity 
potential on the_i_. . ... - 
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The Aerodynamic Forces on the Tail Surfaces 

D erivat ion of t he span load distribution . - As has 
"been mentioned, the I of figure 2 represents the trail- 
ing vortex sheet of a tail ^roup. The potential function 
just discussed corresponds to the induced flow about the 
vortex sheet when the induced velocity- of the sheet itself 
is unity. The follo-wine; calculations will be nade for 
this case of unit induced velocity behind a tail of dimen- 
sions H, D; the final results, however, will be s;iven 
for any induced velocity w and any tail dimensions h, d 

The circulation per unit length (spanwise) of the 
trailing vortex sheet at a point (x, 0) is the differ- 
ence between the velocities on the upper .and the_ lower 
surfaces at that point, i.e., 

where the subscripts i and a refer to the upper and 
the lower surfaces of the sheet, respectively. The bound 
circulation V at the corresponding point on the vertical 
tail equals the total circulation in the trailing vortex 
sheet between this point (x, 0) and the tip (H, 6), 
i.e.,- 

x x 

or T(x) = tf> 2 (x) - « x (x) (18) 



since v x (E) = § 2 (H) 

The bound circulation at a point (0,_jr)_ on the horizon- 
tal t?>,il is ^iven by a corresponding expression. Equa- 
tions (18), (17), and (16) suffice to determine the span- 
wise distribution of circulation on the vertical tail. 

Load distributions have been calculated for various 
values of r and are given in figures 12 and 13 . ¥h.o 
curves for the vertical-tail surface have zero slope at 
x =-0, since x = 0 (cos 9 = m) is a stagnation point 



H . ' IT 
x 
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of the -ootential flow where:- (1^) * (1^). = 0 and, there- 

\9x/ x ^odc's 

fore , Y = —■ = 0. 

&x 

It will "be observed that figures 12 and 13 have non- 
dimensional coordinates and apply to tails of arbitrary 
dimensions and arbi t rat?' induced' wake velocity. 

The t otal forces .- The total lift, forces on the hori- 
zontal and the vertical surfaces are proportional to the 
integrals of T .across the respective spans.' For the , 
vertical tail surface, - 



I ' T(x) dx = / (Q x - Q a ) dx ._'«= 2 /° 3> x dx 



•J 

o 

61 • 



where 81 = cos" 1 (-m), Qs = 0 , • and the integration is 

along the arc of the unit circle from 8 X to 9- a . The . 

velocity potential $3. is "?iven by equation (17) and 

dx. dj5^ is obtained from equation (16). 
d8 ^ d9 ^ 

There results finally 

n 

I r dx = ft (1 + m) (3 - m) (19-) 



Similarly, for the integral of- T on one-half the hori- 
zontal tai"l— surf ace , ' -' 



/ t. * 1 ,„ v/ v , 3 - m + 2 J 2 (l - m) 

/ Tdy = % (3 - m) (l + m) log 1 

•y <= 1 + m 



1 + m 

- (1 +. m)' y 2 (1 - m) (20) 
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The moment, about the point of intersection, of the 
lift on the vertical tail is proportional to 

Fx ax = 2 / W (z) z dz __- ■ 
o o 

(21) 



3/a r g--i 

(1 ~ m a ) + 4 (l + m) cos" (~m) + a v 1 - m 



Similarly, for the moment of the lift on the horizontal 
tail, ' 

D . 

2 Tydy = (l-m a ) 3/a + 4(l+o) fcos" 1 (-m) + mVl-m 2 J 
o 

- 4tt (l+iii) (22) 

In expressions (21) and (22), cos"" 3 " (-n) is taken in the 
first or the second quadrant. The sum of these expres- 
sions corresponds to the resultant moment, 

H D 
^ Tx dx + 2 / T ydy = 

■"0 0 • - 

(l-m s ) 3/3 + 8(l+m) [cos-^-a) + m </l-m 2 J -4tt(1+ii) (2ET) 

The foregoing relations hold for the case of unit 
normal induced velocity and tail dimensions H and D. 
General formulas, applicable to tails of arbitrary dimen- 
sions h and d and arbitrary induced-walce velocity w, 
are obtained from these relations by multiplying by the 
proper factors and substituting from equations (15). '■_ 
Thus from equation (19), . — 

L h = p V w Qj) ir(l + m) ( 3 - m ) 
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or 
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L h . £ pA* (*•) J (3 - m) 



Similarly, 



IT - COS 



1 m + ffi-s/l « m 



Zj2 (1 + m) 



(a) 



(o) 



U = ipv a (lV j"^io-) 3/3 - nivCT-co a -^ m |f -l__ J! __'] . , 



She ratio w/V that appears in these expressions will oe 
discussed in the following section. 
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The limiting forma of these equations for the forces 
and the moments as the ratio of vertical tail span to hor- 
izontal tail span approaches zero or infinity are of in- 
terest. Thus, for finite h, 

11. Ifc -| PVV « (I) 11. L h - |>h' | (fX - 



r— > 0 



lim = co lim 

r — »-0 r — -oa 



o p — p. oo 



lim M d = o= lim M d = 0 

r — * 0 , r — *► oo 

where 



a o a o 



— f TTA 0 g + g 

In the preceding -expressions for w/V,_ A Q and A^, 

correspond to semi ellipti cal and elliptical chord' distri- 
butions, respectively, as vrill he shown later. 

The normal induced velocity and effective aspect 
r atio . - Inasmuch as w/2V is the induced angle at the 
vertical-tail surface, ' ' " ~ ' , 



O ih - a 0 (-, - Jt) ■ '■ XM) 



Division of equations (24a) T?y ^ PV 2 S h ?ives 



E 



liminating from equations (25). and (26) gives 
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T 



a„ a A 
o o 



S° + - (3 - in) 



4 S v 



(27X 



Substituting for w/2V in equation (25) gives the slope 
of the lift curve: 



= a. 



1 - 



o 



a„ + - — (3 - n) 



(28) 



Now, for a plain elliptical airfoil . (reference 6), 

a = a„ (l -S— -) (29") 

Comparison of equations (28) and (29) gives a value, for 
the effective aspect ratio pf the vertical tail surfaco : 

.2 



or 



A e = -f- (3 - ») 



3 - m 



(30) 



Talues of m 
figure 14. For 



and A /A have "been plotted, against r in 



05 (no end plate), 



♦ 1, the 



effective aspect ratio is the actual aspect ratio, and the 
vertical tail -load distribution is elliptical, For r — -0 
(infinite end plate), ra — -1 , the effective aspect 
ratio is twice the actual aspect ratio, and the vertical- 
tail- load distribution is semielliptical. 

Chord distributing-" The chord distribution: of th.c 
idealized tail-surface combination is related to the do- 
rived span load distribution by the equation 



cVa 



2 (a a - a,) 



Because both a Q and are constant along the span 
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of the vertical .tail , its chord distribution must he pro- 
portional to its load distribution. Thus, the curves of 
figure 12 indicate not only the load "but also the chord 
distributions. It will be seen from the figure that, for 
values of r lying in the useful working range r < 1, 
the vertical-tail chord distributions appropriate closely 
to the seni elliptical. This fact justifies the compari- 
son of results obtained by "this method with those of part 
I (in which the vertical-tail chord distributions were 
chosen to be somielliptical) . ' For the .horizontal tail, 
both oc^ and oc,j_ are -zero; hence its chord nust be infi- 
nite. 

COMPARISON OP TH3 TWO METHODS 



A comparison of results obtained by the two methods - 
is given in figures 5 and 15. It appears that"T;ho implied 
infinite-chord horizontal tail for the condition of mini- 
mum induced drag introduces an excessive end-plate ofrfect. 
Thus for the tail group of figure 15, the indicated values 

of Oj. , g l > and ^- e /^- are about 10 percent, 100 per- 
il d 

cent, and 30 percent higher, respectively, than tho values 
obtained by tho lifting-lino-equation method. 

.That the usual method of calculating end-plate effects 
is considerably in error in the present instance is evi- 
dent. It may be noted, however, in figure 5, that the rel- 
ative error decreases with increasing r, that is, with 
decreasing span of the end plate. The matter is of inter- 
est with regard to the theoretical studies of twin tails 
(references 1 and 2). The relative span of the end plates 
(the vertical tails in this case) is so smail that the er- 
ror incurred in using the minimum induced-drag method to' 
determine the increased of fecti vonoss of the horizontal 
tail may well be negligible. 



Lan^ley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, 7a., December 30, 1940. 
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Tlgura 5.- TaarijrtLoii of with r for wftoaa 
tt.aaJ.-tal! aapeet ratloa. 
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[ Flfurs 6.- T)u amtti of tha bortiontal and rartloal tail forota 
about tha tntaraaotioa for Yarioa* Yirtloai-tail 
aapaot ratio*. 
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Figs. 14,15 
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Figure 14.- Variation of A Q /A with r. Minimum 
induced drag. 
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Figure 15.- Comparison of span load distributions 

calculated by the two methods. Vertical- 
tail aspect ratio, 1.375. ; r , 0.33. 



